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We investigate the Q 2 dependence of inclusive electron- 
proton scattering F2 structure function data in both the nu- 
cleon resonance region and the deep inelastic region, at mo- 
mentum transfers below 5 (GeV/c) 2 . From these data we 
construct moments of F2, down to momentum transfers of 
Q 2 « 0.1 (GeV/c) 2 . The second moment is only slowly vary- 
ing with Q 2 down to Q 2 ~ 1 (GeV/c) 2 , which is a reflection 
of duality. Below Q 2 of 1 (GeV/c) 2 , the Q 2 dependence of 
the moments is predominantly governed by the elastic contri- 
bution, whereas the inelastic channels still seem governed by 
local duality. 



I. INTRODUCTION 

Perturbative QCD (pQCD) describes the scale depen- 
dence of a wide number of hard processes, up to Next 
to Leading Order (NLO), and in a few cases up to Next 
to Next to Leading Order (NNLO) [Q. However, direct 
comparisons of calculations with experimental cross sec- 
tions are often affected by the uncertainty coming from 
non-perturbative contributions, that are less well under- 
stood, although exceptional cases exist where the non- 
perturbative corrections can be deduced from the large 
distance behavior of the perturbative series || . To this 
date, the mechanisms by which a pQCD description of 
deep inelastic observables starts failing, giving way to 
non-perturbative (NP) behavior, are still largely unde- 
termined. New studies of kinematic regions, and possible 
observables expected to be most sensitive to this transi- 
tion, are currently being pursued. 

This paper is dedicated to interpreting a somewhat 
surprising behavior of the nucleon structure function F% 
in the nucleon resonance region, as recently measured at 
Jefferson Lab and reported in for four- momentum 

transfers squared 0.45 < Q 2 < 3.3 (GeV/c) 2 . This behav- 
ior is consistent with older, lower precision, SLAC data, 
down to Q 2 = 0.07 (GeV/c) 2 [§. At the lower values of 
Q 2 , the resonance region, typically defined to comprise 
the region in invariant mass M 2 < W 2 (= M 2 + Q 2 (l/x - 
1)) < 4 GeV 2 , with M the proton mass and x the Bjorken 
scaling variable, is at low values of x, x « 0.1. Thus, the 
resonances do not characterize the large x behavior that 
was extensively studied before, in e.g. Refs. ||-^]. It 
was, however, found that the resonance region seems to 



average to a single curve even at these very low values 
of Q 2 . This scaling curve mimics a valence-like quark 
distribution, the new data of || adding new information 
for x < 0.25. 

It is the aim of this paper to present a low-Q 2 mo- 
ment analysis of the inclusive electron-proton scattering 
i<2 structure function data, based upon the new precise 
resonance region and existing deep inelastic region data, 
and to interpret these moments in the context of the sur- 
prising resonance region behavior of |}| . Such an analysis 
would extend previous pQCD analyses of the type per- 
formed in |Io|Jrj-P] , to the lowest values of Q 2 and x of || . 
In particular, as Q 2 decreases, we expect NP effects to 
dominate the cross section, reducing the agreement with 
pQCD predictions. Such an analysis may be complemen- 
tary to the ones addressing another set of low-Q 2 mea- 
surements of F2 obtained recently at HERA p| , which 
however belong to a completely different kinematical re- 
gion, characterized by very low x (down to x ss 10~ 6 ) 
and, accordingly, to very large W 2 . In |ll[] it has been 
found that the structure function Fi can be described 
by pQCD down to Q 2 » 1 (GeV/c) 2 , provided a gluon 
distribution vanishing at low x (also referred to in the 
literature as "valence- like" [^l) ) and a non vanishing but 
small sea distribution, are adopted. However, in Ref. 
pr[ it was also found that, at even lower values of Q 2 , 
0.045 < Q 2 < 0.65 (GeV/c) 2 , and in the same very low 
x region, the pQCD description breaks down giving way 
to different types (e.g. Regge type) of descriptions. 

The fact that pQCD can be extended to very low Q 2 
(Q 2 « 1 (GeV/c) 2 ) was unexpected and supports the idea 
that the QCD coupling constant might have a smoother 
behavior at low Q 2 than predicted by asymptotic free- 
dom expressions. Or, in other words, the confinement 
mechanism might manifest itself in a softer way [(H^-fbil . 
It is also unexpected from the point of view of a pQCD 
based analysis that the resonance data, lying on the low 
side of the invariant mass spectrum, would follow a curve 
that is not too far from the DIS valence quark curve down 
to very low Q 2 . Is this still a signature of pQCD? 

In Section II we remind the reader of the Bloom- 
Gilman duality phenomenon. In Section III we review 
the world inclusive (e,e') data at low Q 2 , both in the low 
and in the large W 2 regions. In Section IV, we evaluate 
the moments of the proton structure function in the low 
Q 2 region, including the elastic contribution, and we ex- 
tend the pQCD based analysis of |j| to this region. In 
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Section V, we discuss the role played by the elastic contri- 
bution in rendering the Q 2 dependence of the moments. 
In Section VI, we draw our conclusions. 

II. BLOOM-GILMAN DUALITY 

Three decades ago, Bloom and Gilman observed a fas- 
cinating correspondence between the resonance electro- 
production and deep inelastic kinematic regimes of inclu- 
sive electron- nucleon scattering |D| . Specifically, it was 
observed that the resonance strength could be related to 
the deep inelastic strength via a scaling variable which 
allowed for a comparison of the lower W 2 and Q 2 reso- 
nance region data to the higher W 2 and Q 2 deep inelastic 
data. It was observed that the deep inelastic data are ap- 
parently equivalent to an average of the resonance region 
data. Furthermore, this behavior was observed over a 
range in Q 2 and W 2 , and it was found that, with chang- 
ing Q 2 , the resonances move along, but always average to, 
the smooth scaling curve typically associated with deep 
inelastic scattering (DIS). This behavior clearly hinted 
at a common origin for resonance (hadron) electropro- 
duction and deep inelastic (partonic) scattering, termed 
parton-hadron, or Bloom-Gilman, duality. 

A global kind of parton-hadron duality is well estab- 
lished: low-energy resonance production can be shown to 
be related to the high-energy behavior of hadron-hadron 
scattering ]l6|-p^[; the familiar ratio of e + e~ — > hadrons 
over e + e~ — > muons uses duality to relate the hadron pro- 
duction to the sum of the squared charges of the quarks: 
here duality is guaranteed by unitarity (in this, one could 
argue that the p production channel exhibits local dual- 
ity, in that its area averages to about the same global 
value) [ fL9f ; in pQCD the high- momentum transfer be- 
havior of nucleon resonances can be related to the high- 
energy transfer behavior of DIS However, it is 
not clear why duality should also work in a localized re- 
gion, and even at relatively low momentum transfers. 

Inclusive deep inelastic scattering on nucleons is a 
firmly-established tool for the investigation of the quark- 
parton model. At large enough values of W 2 and Q 2 , 
with W 2 >> Q 2 , a precise description of the Q 2 be- 
havior of the nucleon structure function Fi = VW2 can 
be given in terms of a perturbative series in as (Q 2 ), up 
to NLO pi] , p^| . Such Q 2 behavior becomes especially 
transparent in comparing the high Q 2 (> 10 (GeV/c) 2 ) 
moments of F2 with pQCD predictions poUl- 

An analysis of the resonance region at smaller W 2 and 
Q 2 in terms of QCD was first presented in Refs. p0[ , 
where Bloom and Gilman's duality was re-interpreted: 
The integrals of the structure function, performed in |lq ] 
over the energy transfer is, were translated into inte- 
grals over the variable x (or Nachtmann £ = 2x/(l + 
y/l + AM 2 x 2 /Q 2 p3[ ), in order to account for finite tar- 
get mass effects). Bloom-Gilman duality was translated 



into a correspondence between the n —2 moment of the 
i*2 structure function in the low Q 2 region, characterized 
by resonances, and in the high Q 2 scaling region, respec- 
tively. The fall of the resonances along a smooth scaling 
curve with increasing Q 2 was to be attributed [jl0| to 
the fact that there exist only small changes in the low n 
moments of the F2 structure function due to power cor- 
rections in addition to the predicted perturbative ones. 
The appearence of power corrections is interpreted as a 
signal of deviations of the inclusive cross section from 
perturbative predictions, which one can envisage as due 
to the increasing importance of interactions between the 
quark struck in the electron-nucleon hard scattering pro- 
cess and the other quarks in the nucleon. In inclusive 
DIS, the Operator Product Expansion applies and power 
corrections are determined by the matrix elements of op- 
erators of higher twist (dimension-spin), which are re- 
lated to multiparton configurations. 

Such effects are inversely proportional to Q 2 , and can 
therefore be large at small Q 2 . If they are not, averages 
of the F 2 structure function over a sufficient range in x at 
moderate Q 2 are approximately the same as at high Q 2 . 
Notwithstanding, the dynamical origin of local duality, 
and thus the reason why the higher- twist contributions, 
undoubtedly required to construct the coherent nucleon 
resonances, tend to largely cancel on average, even at 
momentum transfers below 5 (GeV/c) 2 , is still not un- 
derstood ||||. 

III. LOCAL DUALITY AT LOW MOMENTUM 
TRANSFER 

In Fig. 1 we show an overview of recent high-precision 
proton resonance F2 data at low Q 2 pj. We also in- 
clude data from SLAC at Q 2 < 0.3 (GeV/c) 2 f§. For 
the former, the systematic uncertainty is estimated to 
be 3.5% [Q. For the latter, due to uncertainties in abso- 
lute normalization and radiative corrections, we estimate 
the systematic uncertainty to be better than 10%. The 
solid curves represent, for the different kinematics, the 
single scaling curve defined by averaging all nucleon res- 
onance F 2 data, regardless of Q 2 , W 2 , as a function of £ 
0. As one can see the individual spectra, at various Q 2 , 
oscillate around this single-curve parameterization. We 
emphasize that this is not by construction, as the pa- 
rameterization, at any given value of £, is obtained from 
a range of nucleon resonance data at variant values of Q 2 
and W 2 (e.g., the second resonance bump could have al- 
ways been below the scaling curve, while the first above, 
etc.). Our main observation is that apparently nature 
forces the oscillatory behavior of the various resonance 
bumps around a scaling curve, which has a valence-like 
behavior. 

This has been studied quantitatively in Q where it 
was observed that the behavior of averaged nucleon res- 
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onance data at £ > 0.3, corresponding to Q 2 > 0.5 
(GeV/c) 2 in the nucleon resonance region, is indistin- 
guishable from the F% DIS behavior, consistent with the 
findings of Bloom and Gilman |15|. The behavior of av- 
eraged nucleon resonance data for £ < 0.3, corresponding 
to Q 2 < 0.5 (GeV/c) 2 in the nucleon resonance region, 
mimics Q xF% data obtained from averaging neutrino 
and antineutrino DIS data [Q. The latter, to leading or- 
der in QCD, selects the difference of quark and antiquark 
distribution functions, and is predominantly sensitive to 
a valence quark only distribution. 

Increasing from Q 2 w 0.07 (GeV/c) 2 (Fig. la) to Q 2 w 
3 (GeV/c) 2 (Figs. lh,li) the F 2 spectra change shape 
drastically. The low Q 2 spectrum shows a predominant 
N—A transition (we do not show the elastic peak, huge at 
this Q 2 ), and relatively minor strength at larger energy 
transfers. This is not surprising, as at these relatively 
small energy and four-momentum transfers one would 
expect to predominantly excite the valence quarks. At 
Q 2 = 3.0 (GeV/c) 2 , on the other hand, one sees that 
the prominent nucleon resonances are largely reduced, 
and the inelastic background enhanced. Furthermore, F2 
in the higher resonance regions is larger than F2 in the 
N — A transition region. Apparently, a swap of strength 
has occurred between the various channels as a function 
of Q 2 . 

To further illustrate how the nucleon resonances seem 
to follow a valence-like curve, we show in Fig. 2 the be- 
havior of the N — A transition region (here defined as 
1.2 < W 2 < 1.9 GcV 2 ) and the second resonance region 
(defined as 1.9 < W 2 < 2.5 GeV 2 ) for various Q 2 as a 
function of £, in comparison with the global scaling curve 
defined in |1, and used in Fig. 1. As concluded in Refs. 
HQ , it seems that the nucleon resonances slide along one 
global scaling curve (note that the apparent difference in 
scaling curve between Figs. 1 and 2 reflects only the con- 
version from W 2 to £, for fixed Q 2 )- One can see that, 
if nature forces the oscillatory behavior around a global 
scaling curve even at low Q 2 , the resonance excitation 
strengths will necessarily grow in the region below £ ss 
0.25 where the maximum of the global scaling curve oc- 
curs, and subsequently decrease once the maximum of 
the global scaling curve has been crossed. Compare, for 
instance, the behavior of the N — A transition region with 
that of the larger-mass resonance regions: at Q 2 = 0.45 
(GeV/c) 2 (solid circles in Fig. 2 (top), and Fig. lc) the 
N — A transition region strength is large being at about 
the maximum of the scaling curve. Its strength, as for all 
Q 2 < 0.45 (GeV/c) 2 , is also larger than the higher-mass 
resonance regions which lie at lower £ for the same Q 2 . 
On the other hand, for Q 2 rj 3.0 (GeV/c) 2 (open circles 
in Fig. 2 (top), and Fig. lh) the N — A transition region 
strength is small because it is positioned at large £, and 
smaller than the higher-mass resonance regions that lie 
at lower values of £ , but have crossed the maximum of the 
scaling curve. The smooth curve, to which the nucleon 



resonance regions tend, determines the momentum trans- 
fer dependence of the various nucleon resonance regions, 
forcing the nucleon elastic and transition form factors to 
scale like Q~ 4 |2f| at relatively small Q 2 , resembling the 
Q~ 4 scaling as predicted by QCD counting rules. 

We now ask the question: Can the observed behavior 
of the averaged nucleon resonance spectra be explained 
within pQCD? 

In Fig. 3, we show a compilation of the world's 
electron-proton scattering data for the F2 structure func- 
tion at low Q 2 . The deep inelastic (W 2 > 4 GeV 2 ) data 
originate from SLAC @), CERN (NMG) |27), FNAL 
(E665) @, and DESY (HI, ZEUS) @-|||l|! As be- 
fore we include data in the proton resonance region from 
SLAC || and JLab [pi. The solid curves indicate the 
next-to-leading order parameterizations of Gliick, Reya, 
and Vogt (GRV) |34|,|35|], which use input parton distri- 
butions starting from very low Q 2 values. In the third 
panel down of Fig. 3, the GRV calculation is the GRV in- 
put distribution at Q 2 = 0.4 (GeV/c) 2 J35|, without any 
evolution, neglecting sea and gluon contributions. As 
one can see, the proton resonance region data for Q 2 = 
3.1 (GeV/c) 2 and Q 2 = 0.9 (GeV/c) 2 smoothly join the 
deep inelastic data, and agree well with the GRV next- 
to-leading order calculations, exhibiting the local duality 
witnessed by Bloom and Gilman |l5|| . 

Turning our attention to the bottom two panels of Fig. 
3, the low Q 2 F2 data, we are only left with the recent 
DESY data [§2],||, some sparse FNAL data ||, and the 
nucleon resonance data §1). The DESY data exhibit the 
well-known collapse of the proton structure functions at 
(very) small x. One interpretation of this effect is that, 
at these small Q 2 , one sees a smooth transition from DIS 
to the real photon point at Q 2 
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Gauge 
= 0, the 



|3| 

invariance requires that, for consistency near Q 2 
structure function F 2 for inelastic channels must vanish 
like Q 2 a( 7 p)/(47r 2 a em ) ||. 

We emphasize here the difference in reaching the low 
Q 2 region for the various values of x. In Fig. 4 we dis- 
tinguish among three different limits (indicated by the 
three arrows in the figure). For the DESY experiments, 
low Q 2 is established at small x by having a large (« con- 
stant) amount of energy transfer v. For this reason, it is 
expected that this region exhibits similar characteristics 
as the parton model. For the JLab/SLAC experiments 
at x w 0.1, one reaches low Q 2 at relatively small en- 
ergy transfers. If the limit is taken at fixed x, x « 0.1, 
the resonance data do not exhibit such a drastic collapse 
with Q 2 as observed at HERA, and they stay fairly con- 
stant. In fact, they still seem to oscillate around one 
Q 2 -independent global curve p], informing us that the 
Q 2 dependence of the larger x nucleon resonance data is 
rather shallow. It is this behavior that makes the scal- 
ing curve of Figs. [1,2], defined by the world's nucleon 
resonance F2 data, look "valence-like" ||. If, instead, 
one takes the limit at a fixed W 2 , by restricting oneself 
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always to the nucleoli resonance region (M 2 < W 2 < 
4 GeV 2 ), then one can see that the nucleon resonances 
"slide" down the x scale to lower x for lower Q 2 (e.g. see 
the constant W — 2 GeV arrow in Fig. 4), where their 
strength dies out as a function of Q 2 . 



IV. 



MOMENTS OF F 2 P 



We construct the experimental moments of the struc- 
ture function, F 2l for the Q 2 range up to 10 (GeV/c) 2 . 
The Cornwall- Norton moments are defined as 
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M n (Q z )= / dxx n - 2 F 2 (x,Q'), 
Jo 

and the Nachtmann moments as 
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3 + 3(n + l)r + n(n + 2)r 2 



(n + 2)(n + 3) 



vW 2 {x,Q 2 ). (2) 



Here, r = (1 + AM 2 x 2 /Q 2 ) 1 ^ 2 , and x t hr is Bjorken x 
for pion threshold. We add to these integrals the elastic 
contribution, at x = 1, where 



G 2 E (Q 2 ) + ^G 2 M (Q 2 ) 
vW 2 {x,Q 2 ) =(5(1 -x)± -. , (3) 
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Ge (Gm) is the proton electric (magnetic) form factor. 
For the proton form factors, we use a fit to the world's 
data by Bosted Q. 

To obtain the inelastic contributions, we integrate data 
like those shown in Figs. 1 and 3. Apart from the data 
shown in Fig. 1, we have added Q 2 points where addi- 
tional data were available. For Q 2 < 0.6 (GeV/c) 2 , we 
have constrained our search to elastic and nucleon reso- 
nance data. For 0.6 < Q 2 < 4 (GeV/c) 2 , we have used 
nucleon resonance data in combination with deep inelas- 
tic data, whereas for Q 2 > 4 (GeV/c) 2 we have con- 
structed the moments utilizing both deep inelastic and 
nucleon resonance models, similar as in Ref. ||. For the 
smallest values of Q 2 (< 0.6 (GeV/c) 2 ), we assume a con- 
stant value of F 2 below x for W 2 — 4.0 GeV 2 , as only 
sparse data exists. As one can see from Fig. 3, this may 
not be a bad approximation for Q 2 < 0.6 (GeV/c) 2 , es- 
pecially since the nucleon resonance region data extend 
down to x < 0.1, and the integration area below x = 0.1 
is expected to be small only. To judge the uncertainty 
in this procedure, we have also integrated the Q 2 « 0.2 
(GeV/c) 2 data starting at W 2 = 9.0 GeV 2 (rather than 
W 2 = 4.0 GeV 2 ). This changes the second moment by 
less than 3%. Lastly, in some cases, we used a model to 
construct data at fixed Q 2 , rather than allowing for the 



small range of Q 2 in the data. This effect on the second 
moments was found to be small, < 3%, and far smaller 
for the higher moments. Thus, we believe the total uncer- 
tainty in the moments we calculate to be less than 5%. 
We show the values for the second, fourth, sixth, and 
eigth Cornwall-Norton (top) and Nachtmann (bottom) 
moments of the proton, extracted from the world's data, 
including deep inelastic, nucleon resonance, and elastic 
data, as described above, in Fig. 5. Similarly, Tables I 
and II list the numerical values of the moments, with the 
elastic contribution to each separately given. 

As expected, the elastic contribution dominates the 
moments at the lowest Q 2 . Note that the Cornwall- 
Norton moments will become unity, i.e. the proton 
charge squared, at Q 2 = 0, whereas the Nachtmann mo- 
ments will vanish at Q 2 — 0, as can readily be seen from 
Eqn. 3. This can be attributed to the fact that, with 
respect to Bjorken x, the Nachtmann scaling variable £ 
correctly takes into account the finite proton mass scale 
|p3| , but does not account for any other significant mass 
scale (like the quark masses). As the interpretation of 
the Cornwall-Norton moments in the Q 2 < 1 (GeV/c) 2 
region seems more intuitive, and we are interested here in 
the low-Q 2 behavior of the moments, we will concentrate 
on these moments in the remainder of the discussions in 
this work. To emphasize that there is indeed not much 
difference between the Cornwall-Norton and Nachtmann 
moments if one concentrates on the low-Q 2 region where 
the elastic contribution turns dominant, Fig. 6 graphi- 
cally displays the relative contribution of the elastic chan- 
nel to the total moment for both Cornwall-Norton (top) 
and Nachtmann (bottom) moments, for n = 2 (solid cir- 
cles), n — 4 (squares), n = 6 (triangles) and n = 8 
(stars), from Tables I and II. Nonetheless, as the ben- 
efit of the Nachtmann moments is to push an analysis in 
terms of an Operator-Product Expansion to lower values 
of Q 2 , taking correctly target-mass effects into account, 
we show everywhere similar figures, for comparison, for 
the Nachtmann moments. Note that one can argue that 
the relative contribution of the elastic grows slower for 
Q 2 — > if one uses the Nachtmann moments. 

First, let us revisit the large Q 2 , Q 2 > 10 (GeV/c) 2 , be- 
havior of the moments. At asymptotically large Q 2 , one 
can write the non-singlet moments M^ S (Q 2 ) at Leading 
Order (LO) [] in the perturbative expansion, as p|Lp| 



M^ s (Q 2 )=A n (ln(Q 2 /A 2 ))- 1 /<* 
where A is the QCD scale parameter and 



d» s = 7^/2/3, 



(4) 



(5) 



*NLO corrections are important in general and will be dis- 
cussed in detail in p§]. Their discussion is, however, not 
essential for discussing the points raised in the present paper. 
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where f3 a = 11 - (2/3)Nf, with Nf the number of flavors, 
and 7o, n Ar,s the leading-order non-singlet anomalous di- 
mensions numerically specified in p2fl . To circumvent 
the requirement of non-singlet moments we highlight in 
Fig. 7 the n = 4 moment: although we use the same F 2 
structure function data for all moments, the weighting 



with 



in these moments will emphasize the large- a; 



region, at higher n, and thus approximate a non-singlet 
moment. On the top of Fig. 7 we show the Cornwall- 
Norton moment, on the bottom the Nachtmann moment 
(each to the power — l/d^ s ). We show the data in a log- 
log plot, in order to encompass also the low-Q 2 region. 
The moments are shown both with (stars) and without 
(open circles) the elastic contribution included. 

The dashed curves exhibit a fit to the data, from 0, 
limited to Q 2 > 20 (GeV/c) 2 to minimize the effect of 
higher order corrections, in the form [M^ fS Y~ 1 ^ dn ) w 
Pi^(Q 2 /A 2 ). In §, this fit gives Pi = 27.46 (27.05) ± 
0.25 (0.24) and A = 250 MeV, for the Cornwall-Norton 
(Nachtmann) moment, rendering the expected logarith- 
mic scaling behavior in QCD at asymptotic Q 2 . The dot- 
dashed (dotted) curves exhibit a similar fit in the form 
(Pi+ P2/Q 2 + P 3 /Q 4 )- 1/d "ln(Q 2 /A 2 ) down to Q 2 = 2.0 
(GeV/c) 2 from the same Reference 0, thus taking into 
account power corrections of order 1/Q 2 and 1/Q 4 (1/Q 2 
only). Numerical values for the 1/Q 2 and 1/Q 4 coeffi- 
cients are P 2 = 0.33 ± 0.04 (0.33 ± 0.04) and P 3 = 4.69 
± 0.19 (1.61 ± 0.15) for the Cornwall-Norton (Nacht- 
mann) moment (see also the caption of Fig. 7) . One can 
easily verify from Fig. 7 that the magnitude of the P3 co- 
efficient is in this case dominated by the inclusion of the 
clastic contribution. Similarly, the Nachtmann (n = 4) 
moment analysis gives a drastically different value for 
P3 from the Cornwall-Norton (n — 4) moment analysis 
mainly due to the different contribution from the elastic. 

In order to illustrate the relative strength of the elas- 
tic contribution compared to other W 2 regions, we show 
in Fig. 8 the second (Fig. 8a), fourth (8b), sixth (8c) 
and eigth (8d) Cornwall-Norton moments for Q 2 < 5 
(GeV/c) 2 , separated in the elastic contribution (squares, 
due to our choice of vertical scale sometimes only visible 
at the higher Q 2 ), the contribution of the N — A tran- 
sition region (triangles, 1.2 < W 2 < 1.9 GeV 2 ), of the 
second resonance region (open circles, 1.9 < W 2 < 2.5 
GeV 2 ) and of the "deep inelastic" region (stars, W 2 > 4 
GeV 2 ). The total moment is given by the solid circles, 
and the curves connect the various data to guide the 
eye. The chosen finite W 2 regions will start contributing 
to the moments at low Q 2 , recovering part of the loss of 
strength due to the fall-off of the elastic contribution, and 
then also die off, as the resonances move to the larger £ 
side of the scaling curve. The contribution of the W 2 > 
4 GeV 2 region does not die off, as this is not a finite 
W 2 region, so higher- W 2 resonances and/or higher- W 2 
inelastic background start becoming important with in- 



creasing Q 2 , eventually yielding the logarithmic behavior 
of the moments prescribed by QCD. As evidenced by the 
difference between the W 2 > 4 GeV 2 contribution and 
the total moment, the contribution of the region of W 2 < 
4 GeV 2 is non-negligible up to Q 2 m 5 (GeV/c) 2 , even for 
the second moment. Similar remarks hold for the various 
Nachtmann moments, shown in Fig. 9. 

We note here also that the behavior of the second 
Cornwall-Norton F2 moment we extract is very similar to 
the behavior in the second moment of the spin-dependent 
gi structure function ]39[ | . Presently we only have sparse 
<7i data in the nucleon resonance region, for 0.1 < Q 2 < 5 
(GeV/c) 2 , such that we can not verify precisely whether 
the spin-dependent nucleon resonance data tend to oscil- 
late around a similar smooth curve. However, the limited 
data are not inconsistent with such a behavior Eo). Also, 
we presently do not have enough data for the longitudinal 
structure function Fl to verify a similar onset of duality 
iftlfl , although sparse hints do exist in the present world's 
data (42 1 . 



V. DISCUSSION OF RESULTS 

Our findings, that the moments of P2 show a smooth 
transition from DIS down to Q 2 w (GeV/c) 2 , and 
that the nucleon resonances tend to oscillate around one 
smooth curve, support the findings of Ref. p3] , [l4] . p8[ that 
higher-twist effects are small if one looks at the low-Q 2 
behavior of F 2 for Q 2 ~ 1 (GeV/c) 2 . The dynamical pro- 
cess of local duality dictates minimal Q 2 dependence of 
F2 at small Q 2 ; in terms of pQCD this can be explained 
if the higher-twist effects are reduced on average in the 
nucleon resonance region p5fl . Nonetheless, higher- twist 
effects must be responsible for the nucleon resonances 
themselves. The results for the lower moments of P2, 
presented here, show a forced transition from the elastic 
point to the large Q 2 limit, supported by the oscillations 
of the nucleon resonance region around one smooth curve 
at low Q 2 . This smooth curve resembles the deep inelas- 
tic data at Q 2 ~ 1 (GeV/c) 2 , and higher- twist effects 
continue to be small from there on. 

The extension of a pQCD analysis to very low values of 
Q 2 , i.e. Q 2 — > A 2 , is of doubtful predictivity. Here even 
the low n moments are mainly sensitive to the inclusion of 
the elastic contribution (as this contributes already close 
to 10% to the n = 2 moment at Q 2 = 2 (GeV/c) 2 ): a pic- 
ture in terms of non-interacting quark and gluon degrees 
of freedom is clearly no longer tenable. Similarly the most 
plausible interpretation of the behavior observed for the 
dynamical power corrections, which should represent the 
initial signature of NP effects overshadowing the pertur- 
bative expansion, would in fact be that their contribution 
in the 1/Q 2 expansion is enhanced at decreasing values 
of Q 2 . However, this contribution can be damped on av- 
erage if coefficients of different powers are large but with 
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opposite signs, thus causing the cancellation required to 
make the moment a slowly varying function of Q 2 only. 
These cancellations are evidently of NP origin and are 
at low Q 2 not easily associated to partonic degrees of 
freedom in the multiparton correlations. 

One has to be careful making further concrete state- 
ments. E.g., it is often remarked that duality obviously 
does not work universally, as the n — 2 Cornwall-Norton 
moment at Q 2 — > equals the coherent sum of the quark 
charges squared, which can not equal the parton model 
expectation in the Bjorken limit. First of all, F 2 is not 
purely transverse, which casts doubt on drawing too def- 
inite a conclusion at low values of Q 2 . Results of an 
experiment to measure R = <jlI&t in the nucleon res- 
onance region are needed and forthcoming p6| . More 
importantly, the moment analysis will at lower and lower 
Q 2 be sensitive to a smaller and smaller region in £, as 
the inelastic region and nucleon resonances move to lower 
£ and their absolute contribution to the moments dimin- 
ishes. Even in a world without elastic electron-proton 
scattering, one would, at the lowest Q 2 , be mainly sen- 
sitive to a single resonance transition region, in this case 
the N — A transition region. This is illustrated in Fig. 
10, where we repeat Fig. 5 for a world without elastic 
electron-proton scattering. It can easily be seen (Figs. 1, 
3, 8, and 9) that the N — A transition contribution, de- 
fined by the region 1.2 < W 2 < 1.9 GeV 2 , becomes more 
and more dominant at low Q 2 . Thus, in this world a 
moment analysis at low Q 2 would become dominated by 
the coherent N — A process. Still, on average the N — A 
transition region seems to obey duality, in that it oscil- 
lates around a similar smooth scaling curve. This is not 
inconsistent with the findings of |L0]], the discrete res- 
onance transitions are reminiscent of large higher-twist 
effects, but on average still seem to cancel to large extent 
(although quantitatively perhaps not as well as at higher 
Q 2 ). The above example is not dissimilar to stating that 
duality will not work on top of any given resonance peak. 
A region over which to average is always required. 

Furthermore, if one neglects the elastic channel, one 
will at low Q 2 be mainly sensitive to the imposed con- 
straint by gauge invariance that the structure function F 2 
must behave like Q 2 a(jp)/(4ir 2 a em ) |3(|. At the values 
of x where the nucleon resonances are visible at low Q 2 
in Fig. 3 (e.g. x » 0.1), the F 2 structure function does 
not linearly vanish with Q 2 yet, as shown in ||. Thus, 
although the F2 strength in the nucleon resonance region 
has to disappear linearly with Q 2 below some Qq, one can 
argue that the behavior of the data is not reflecting this 
Q 2 < Qq expectation yet. This indicates that the oscil- 
lations the nucleon resonances exhibit around a smooth 
curve, even down to Q 2 ks 0.1 (GeV/c) 2 , is non-trivial. 
As the low-Q 2 F 2 data below W 2 — 4 GeV 2 predomi- 
nantly consists of excited nucleon resonances and hardly 
contributions from inelastic non-resonant processes, one 
can argue that such a smooth curve must be close to a 



curve consisting of valence strength only. In fact, the 
Q 2 dependence of the integrated valence quark strength 
in the GRV model |34||3f| is close to the Q 2 dependence 
of the second Cornwall-Norton moment of F 2 . However, 
this Q 2 dependence is predominantly due to the inclu- 
sion of the elastic channel. Thus, for a picture such as 
the GRV model to be valid, there must be a separate Q 2 
dependence for the vanishing of the large- x strength at 
small Q 2 (governed by the nucleon resonances) and the 
growth of the small- x sea. 

Such ideas are very similar as to what has been ob- 
served in hadron-hadron scattering jlTj. Here, a gener- 
alization of the duality picture was introduced, in which 
resonances were dual to non-diffractive Regge pole ex- 
changes, with the non-resonant contributions dual to 
Pomeron exchange. Within QCD, this corresponds to 
a picture where resonances are dual to valence quarks, 
while the non-resonant background is dual to sea quarks. 
This supports the importance of additional detailed stud- 
ies of electron scattering in the nucleon resonance region, 
in a wide range of x and Q 2 . 

VI. CONCLUSIONS 

We show that the world's data on F 2 , down to Q 2 ~ 
1 (GeV/c) 2 , are reasonably well described within pQCD. 
This includes the nucleon resonance data, which aver- 
age to an approximate scaling curve, due to local dual- 
ity. Down to Q 2 ~ 0.1 (GeV/c) 2 , the nucleon resonance 
data still tend to average to such a curve. The contri- 
bution of the nucleon resonances to the lower moments 
of F 2 dies out at very small Q 2 as they have moved to 
smaller Bjorken x. Instead, the moments below Q 2 rs 1 
(GeV/c) 2 are dominated by the elastic contribution at 
large x. Therefore, a pQCD-based analysis of the low-Q 2 
moments of F 2 will render coefficients of the higher-twist 
terms which are predominantly due to the elastic con- 
tribution. Local duality seems to hold down to at least 
Q 2 w 0.5 (GeV/c) 2 and to prescribe the transition from 
the kinematic region dominated by the elastic contribu- 
tion to the region dominated by deep inelastic scattering. 
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TABLE I. Cornwall-Norton Moments for n = 2, 4, 6, and 8 
at 0.15 < Q 2 < 4.3 (GeV/c) 2 , as extracted from the data (see 
text). The elastic contribution is given as a separate column. 
The uncertainties of the total moments are smaller than 5%. 



Q 2 (GeV/c) 2 


elastic 


n = 2 


n = 4 


n — 6 


n = 8 


0.15 


0.592 


0.652 


0.594 


0.592 


0.592 


0.20 


0.504 


0.584 


0.508 


0.505 


0.504 


0.45 


0.249 


0.379 


0.261 


0.251 


0.250 


0.55 


0.195 


0.341 


0.210 


0.198 


0.196 


0.85 


0.103 


0.278 


0.122 


0.107 


0.104 


0.94 


0.087 


0.264 


0.107 


0.092 


0.088 


1.40 


0.040 


0.231 


0.064 


0.047 


0.043 


1.70 


0.026 


0.219 


0.051 


0.034 


0.029 


2.40 


0.011 


0.203 


0.036 


0.019 


0.014 


3.00 


0.006 


0.196 


0.030 


0.013 


0.009 


3.30 


0.005 


0.192 


0.028 


0.012 


0.008 


4.30 


0.002 


0.184 


0.023 


0.008 


0.005 
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TABLE II. Nachtmann Moments for n — 2, 4, 6, and 8 at 
0.15 < Q 2 < 4.3 (GeV/c) 2 , as extracted from the data (see 
text). The elastic contributions, different for each n, are given 
as a separate entity in the columns. The uncertainties of the 
total moments are smaller than 5%, and the numbers quoted 
can be used to this precision. 



0.15 
0.20 
0.45 
0.55 
0.85 
0.94 
1.40 
1.70 
2.40 
3.00 
3.30 
4.30 



n = 2 
elas. 



total 



0.274 
0.256 
0.160 
0.131 
0.076 
0.065 
0.032 
0.022 
0.010 
0.005 
0.004 
0.002 



0.333 
0.322 
0.281 
0.268 
0.243 
0.235 
0.217 
0.209 
0.197 
0.192 
0.189 
0.181 



n = 4 
elas. 



total 



0.040 
0.047 
0.050 
0.046 
0.034 
0.030 
0.018 
0.013 
0.007 
0.004 
0.003 
0.001 



0.041 
0.049 
0.057 
0.056 
0.047 
0.045 
0.036 
0.033 
0.027 
0.024 
0.023 
0.020 



n — 6 
elas. 



total 



0.0051 
0.0074 
0.0139 
0.0146 
0.0137 
0.0130 
0.0093 
0.0073 
0.0041 
0.0026 
0.0021 
0.0011 



0.0052 
0.0075 
0.0146 
0.0157 
0.0157 
0.0154 
0.0128 
0.0114 
0.0089 
0.0075 
0.0071 
0.0058 



n = 8 
elas. 



total 



0.0006 
0.0011 
0.0037 
0.0045 
0.0054 
0.0054 
0.0047 
0.0040 
0.0026 
0.0018 
0.0015 
0.0008 



0.0006 
0.0011 
0.0038 
0.0046 
0.0057 
0.0058 
0.0055 
0.0051 
0.0040 
0.0033 
0.0031 
0.0024 



0.4 

0.3 

ta 0.2 

0.1 


0.4 

0.3 

0.2 

0.1 


0.4 

0.3 

0.2 

0.1 
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FIG. 1. F2 Spectrum in the nucleon resonance region as a 
function of W 2 , for values of Q 2 = 0.07 (a), 0.20 (b), 0.45 
(c), 0.85 (d), 1.40 (e), 1.70 (f), 2.40 (g), 3.00 (h) and 3.30 (i) 
(GeV/c) 2 . We have superimposed the results from the scaling 
curve from Ref. to illustrate the behavior of the nucleon 
resonance region with increasing Q 2 . 




FIG. 2. F 2 data for the regions 1.2 < W 2 < 1.9 (top) and 
1.9 < W 2 < 2.5 (bottom) GeV 2 , as a function of Nachtmann 
£. Data are shown for Q 2 = 0.07, 0.20, 0.45, 0.85, 1.4, 2.4, and 
3.0 (GeV/c) 2 (left to right), respectively. The solid curve rep- 
resents the scaling curve, determined by averaging all nucleon 
resonance data yfy. 



8 




1.5 



0.5 



0.5 


1.5 



0.5 - 



10 



10 



: Q 2 -3.1 2 - 
- (GeV/c) 


GRV 98, Q; = 3.0 (GeVM 

* JLab 98, Q" = 3.0 (GeV/c)" 

* HI 96 » ZEUS 98 
i ZEUS 95 * E665 96 
' ZEUS 95 o NMC97 
<4* • ZEUS 96 + SLAC 90/92 

1 \ 


: Q 2 -0.9, 
i (GeV/c) 2 

~ A A 

— . 


GRV 98, Q, 2 = 0.9 (GeV/c) 2 

* JLab 98, Q" = 0.85 (GeV/C) 

ZEUS 97 + SLAC 90/92 
» ZEUS 98 
* E665 96 
NMC97 


r q 2 = o.i 

L (GeV/c) 2 
:ai3§l$$ o 


GRV 98, Q 2 = 0.4 (GeV/c) 2 

SLAC, Q 2 = 0.2 (GeV/c) 2 , 
* JLab 98, Q" = 0.45 (GeV/c)" 
ZEUS 97 
■ ZEUS 98 

^ g Q ^ ^ ^^^"^^ \f 


^ Q 2 -o.i 2 

L (GeV/c) 


SLAC, 0.06<Q 2 <Q,09 (GeV/c) 2 , 
ZEUS 97, 0. 1 1<Q"<0. 15 (GeV/c) 




i i «^-*A. 



10 



10 



10 



> 

CD , 















\ 




M^kzr, , ... 







2Mv 



FIG. 4. Kinematics region of the world's F% structure func- 
tion data, for Q 2 < 10 (GeV/c) 2 . Thin lines are for fixed 
Bjorken x. The dashed line is for fixed W = 2 GeV, indi- 
cating the border of the region typically associated with the 
nucleon resonances (the other border being fixed W = M, or 
x — 1). The shaded area exhibits the kinematics region of the 
recent DESY measurements p9|-p3[. The thick arrows indi- 
cate various manners in which one can reach the limit Q 2 — » 
0, at fixed x, fixed W , or fixed v. 



FIG. 3. F% as a function of x for four values of Q 2 , with 
a logarithmic x scale. The symbols indicate various experi- 
ments, as cited in the text. The solid curves in the top two 
panels represent the calculated distributions from the GRV 
collaboration evolved from Q 2 = 0.4 (GeV/c) 2 . The 

solid curve in the third panel represents the input distribution 
itself. 
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FIG. 5. Cornwall-Norton moments (top) and Nachtmann 
moments (bottom) extracted from the world's electron-proton 
scattering data, for n = 2, 4, 6, and 8. The solid lines in- 
dicate the elastic contribution. At low Q (< 4.3 (GeV/c) ) 
the moments (stars) are directly constructed from the world's 
electron-proton Fi database (see text). At larger Q 2 , the mo- 
ments have been extracted from appropriate fits to the world's 
data on inclusive scattering to both the nucleon resonance and 
deep inelastic regions (dashed lines) M. 
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FIG. 7. Log-log plot of the n = 4 Cornwall-Norton (top) 
and Nachtmann (bottom) M n (Q 2 )~ 1 ^ dn moment as a func- 
tion of Q 2 , where d n = ^ n /2f3 pCfl . Stars (open circles) 
do (do not) include the elastic contribution. In the top 
plot, the dashed curve is a fit to the n — 4 moment in 
the form P 1 ln{Q 2 /A 2 ), from |], with Pi = 27.46 (± 0.25), 
and A = 250 MeV. Similarly, the dotted curve uses a form 
(Pi + P 2 /Q 2 )ln(Q 2 /A 2 ), with P 2 = 0.33 (± 0.04), and the 
dot-dashed curve uses a form (Pi+P 2 /Q 2 +P'i/Q i )ln(Q 2 /A 2 ), 
with P3 = 4.69 (± 0.19). In the bottom plot, the fit param- 
eters are Pi = 27.05 (± 0.24), P 2 = 0.33 (± 0.04), and P 3 
1.61 (± 0.15), respectively 
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FIG. 6. Contribution of the elastic channel to the Corn- 
wall-Norton moments (top) and Nachtmann moments (bot- os 
torn) extracted from the world's electron-proton scattering 
data, for n = 2 (solid circles), 4 (squares), 6 (triangles), and 
8 (stars), up to Q 2 = 5 (GeV/c) 2 . The data are from Tables 
I and II. 
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FIG. 8. Second (a), fourth (b), sixth (c) and eigth (d) 
Cornwall-Norton moments. Contributions due to the elastic 
peak (squares), the regions 1.2 < W 2 < 1.9 GeV 2 (triangles), 
1.9 < W 2 < 2.5 GeV 2 (open circles), and W 2 > 4 GeV 2 
(stars) are separately shown, in combination with the total 
moment (solid circles), as a function of the momentum trans- 
fer. Curves connect the various data, and are to guide the eye 
only. 
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FIG. 9. Second (a), fourth (b), sixth (c) and eigth (d) 
Nachtmann moments. Contributions due to the elastic peak 
(squares), the regions 1.2 < W 2 < 1.9 GeV 2 (triangles), 
1.9 < W 2 < 2.5 GeV 2 (open circles), and W 2 > 4 GeV 2 
(stars) are separately shown, in combination with the total 
moment (solid circles), as a function of the momentum trans- 
fer. Curves connect the various data, and are to guide the eye 
only. 
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FIG. 10. Cornwall-Norton moments extracted from the 
world's electron-proton scattering data, for n = 2, 4, 6, and 
8, without the elastic contribution included. The solid 
lines connect the calculated contributions from the N — A 
transition region (1.2 < W 2 < 1.9 GeV 2 ) at various Q 2 . 
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